Vortex core size in interacting cylindrical nanodot arrays 
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The effect of dipolar interactions among cylindrical nanodots, with a vortex-core magnetic config- 
uration, is analyzed by means of analytical calculations. The cylinders are placed in a, N x N square 
array in two configurations - cores oriented parallel to each other and with antiparallel alignment 
between nearest neighbors. Results comprise the variation in the core radius with the number of 
interacting dots, the distance between them and dot height. The dipolar interdot coupling leads to 
a decrease (increase) of the core radius for parallel (antiparallel) arrays. 

PACS numbers: 75.75.-|-a, 75.10.-b 



I. INTRODUCTION 

Regular arrays of magnetic particles produced by 
nanoimprint lithography have attracted strong attention 
during the last decade. Common structures are arrays of 
wires, cylinders, rings, ^^ and tubes, g,!] 

Such structures can be tailored to display different stable 
magnetized states, depending on their geometric details. 
Besides the basic scientific interest in the magnetic prop- 
erties of these systems, there is evidence that they might 
be used in the production of new magnetic devices or as 
media for high density magnetic recording. 9] In partic- 
ular, two-dimensional arrays of magnetic nanoparticles 
have been proposed as candidates for magnctorcsistive 
random access memory (MRAM) devices. [10, 11, 12] 

Recent studies on such structures have been carried 
out with the aim of determining the stable magnetized 
state as a function of the geometry of the particles. 
0) 0, El) 113, dj [Hi In the case of cylindrical particles, 
three idealized characteristic configurations have been 
identified: ferromagnetic with the magnetization paral- 
lel to the basis of the cylinder, ferromagnetic with the 
magnetization parallel to the cylinder axis, and a vor- 
tex state in which most of the magnetic moments lie 
parallel to the basis of the cylinder. The occurrence 
of each of these configurations depends on geometrical 
factors, such as the linear dimensions and their aspect 
ratio T = H/FL with H the height and R the radius of 
the cylinder, [j, 0, [H, Another issue to be consid- 
ered is the interaction between particles. In this case 
the interparticle distance, D, is the important parame- 
ter. [13, til, [H Usually D is large enough to make the 
exchange coupling between particles negligible, making 
the dipolar interaction a fundamental point concerning 
the magnetic state of the system, [ll, [201 

In this paper we focus on cylindrical particles with di- 
mensions such that, in the absence of an external field, 
the magnetic state is a vortex. For this magnetic config- 
uration the dipolar interaction between the dots is due to 
the existence of a core region 21], 22i] in which the mag- 
netic moments have a nonzero component parallel to the 
cylinder axis. In this case it is relevant to understand how 



the core magnetization is affected by the interparticle in- 
teraction. Recently, Porrati et al. [20j investigated an 
array of dots using micromagnetic simulations. In small 
arrays they observed that the dipolar interaction changes 
the size of the magnetic core. However, and because of 
the use of micromagnetic simulations, larger arrays have 
not been investigated. Therefore, analytical calculations 
are very desirable to compare with experiments. With 
this in mind we examine the behavior of the core in ar- 
rays of dots in the vortex-core magnetic phase, in config- 
urations with parallel and antiparallel alignment between 
the cores. We consider analytical calculations based on 
a continuous description of the dots.[23j 



II. SYSTEM & UNITS 

The basic parameters and variables used in our calcula- 
tions are summarized in figure 1 . We consider cylindrical 
dots with radius R and height, or thickness, H, in square 
arrays with N x N dots and center-to-center lattice spac- 
ing D. The distance between any two dots in the array 
is denoted by S. Whenever necessary, cylindrical coordi- 
nates p and (/> are defined on the dot plane normal to the 
cylindrical z-axis. 

The experimental measurement of the vortex core pro- 
file and core radius is not a simple task, and usually it 
is the full core magnetization, /i^, which is measured. 
21, 24] With this in mind, the vortex core is character- 
ized through the calculation of an effective core radius, 
Ceff, defined as the radius of an effective cylinder, uni- 
formly magnetized along its axis, and whose total mag- 
netic moment. 



MoV = MonHClff, 



(1) 



is the same as produced by the z-component of the vortex 
core, as depicted in figure 1. Here V is the dot volume 
and Mq is the saturation magnetization. 

Two types of magnetic ordering within the array are 
examined: all cores parallel (configuration P) and anti- 
parallel nearest-neighbor cores (configuration AP). This 
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choice is not an arbitrary one, because the P configura- 
tion corresponds to the saturated one and the AP con- 
figuration is the ground state of the array. 
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FIG. 1: An illustration of geometrical parameters defining 
each dot and the array. 

All linear dimensions will be considered in units of the 
exchange length L^, defined as = ^2AI [iqMq. The 
dimensionless geometrical parameters are then defined as 
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III. THEORETICAL MODEL 



Large arrays can be studied if we adopt a simplified 
description of the system in which the discrete distribu- 
tion of magnetic moments is replaced with a continuous 
one, defined by a function M{f) such that M{r)5V gives 
the total magnetic moment within the element of vol- 
ume 5V centered at r. Using such a description, the 
internal energy of the array can be written in terms of 
the self-energies, Eg^if, that is, the energies of the iso- 
lated dots, and the interaction contribution, E^nt^ cor- 
responding to the dipolar coupling between the dots. 
Both EsEif and Eint have a magnetostatic term given by 
Ed^p = (a*o/2) J M -VUdV , with U{f) the magneto- 
static potential. Assuming that M{r) varies slowly on 
the scale of the lattice parameter, the exchange term can 
be approximated by E^x = A J ^(Vmi)^(il^, where mi 
is the i-th component of the reduced magnetization with 
respect the saturation value A/q, that is, = Mi/Mo, 
for i — X, y, z.(23| 



A. Vortex-core magnetization 

For the vortex-core configuration we assume that the 
magnetization is independent of z and (j), that is, 



m(r) — mz{p)z + m^{p)(j) 



(3) 



where p is the radial coordinate, z and 4> are unitary 
vectors in cylindrical coordinates, and the normalization 



specifies the core profile, for which we adopt the model 
proposed by Landeros et aZ,[l3l given by 



mM=\^-{p/Bf 



< p<B 



(4) 



and mz{p) — Q \i B < p < R. Here B is a parameter 
related to the core radius and the exponent n is a non- 
negative integer. Alternative expressions for mz{p) have 
been proposed in the literature. [H, [H, 113, HI] Figure 2 
illustrates the calculated magnetization profile (p) for 
a Fe dot {R — 28.2 nm and H — 37.6 nm) using different 
models. The (gray) dash-dotted line corresponds to the 
model by Usov et al 26], the (blue) thin line corresponds 
to the one proposed by Aharoni,|27[ the (black) thick 
line corresponds to our model[14j with n = 4, the dashed 
(red) line represents mz {p) using the model proposed by 
Feldtkeller et aZ[2^, and the dotted (green) hue repre- 
sents the profile using the model presented by HoUinger 

et o/lii. 
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FIG. 2: Calculated magnetization profiles, m^^p), for different 
core models. Dimensions are R = 28.2 nm, H = 37.6 nm and 
the material (iron) es defined by = 3.327 nm. 

Note that the model proposed by A. Aharoni[27l| is 
equivalent to our model (Eq. 2]) using n = \. We can 
observe that our model with n = A agrees well with the 
ones presented in [25] and [28], and has no discontinu- 
ities at mzip = i?) = 0, as in the models proposed by 
Usov et al^ff\ and Aharoni, [27] making the vortex core 
profile easily integrable. The form for mz{p) presented 
above (Equation ^) allows us to obtain the magneti- 
zation profile by minimizing the total energy of the ar- 
ray. The value of B is determined in the minimization 
process, as explained below, and the effective core radius 
{Ceff} can be evaluated equating the z-component of total 
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magnetic moment (Equation ([1])) with the core magnetic 
moment, given by 



fJ-z ^ Mq / mz{p)pdpdzd(j) . 
Jv 



(5) 



Using the proposed model for mz(p), Equation (|4]), we 
obtain pz = MqttHB'^/ (n + 1) and therefore, 



B 



(6) 



B. Total energy calculation 



We write aU the energies in the dimensionless form, 
E = E/pqMqL^. The energies in the P and AP configu- 
rations will be denoted as E'^ and E~ , respectively, and 
can be written as 



E^ = N^E, 



self ^ 



El 



(7) 



where the first term includes the self-energies of the N'^ 
isolated dots, and the second term is the interdot mag- 
netostatic coupling. 



1. Self-energy 

We consider two contributions to the self-energy, Eseif: 
of an isolate dot in the vortex-core configuration: the 
dipolar, S^ipj ^^nd the exchange, E^.^^ contributions. An- 
alytical expressions for the dipolar and exchange ener gies 
have been previously calculated by Landeros et aL[14| 
The dipolar energy results as 
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where 



and 



-)2n-l 



r(n + 3/2)r(2n-F5/2) 



(n + 1)2 



F„(x) =p Fq[(1/2, l,n + 3/2), (n + 2, 2n + 3),-4x^] . 

Here, pFq denotes the generalized hypergeometric func- 
tion, r (x) is the gamma function and n corresponds to 
the exponent used in the core model defined by Equation 

For integer values of n, the exchange energy reduces to 




Eex = tt/i [ln(r/6) +7,, 



(9) 



where 7„ = (l/2)7^(2n) - n7^(-l/2n). Here, H (z) = 
X^i^i [1/* ~ iZI* "I" •^)] generalized harmonic num- 

ber function[29| of the complex variable z. 



2. Interdot magnetostatic coupling 

The dipolar interaction between any two dots in the 
array depends on the center-to-center distance between 
them, 5*, and on the relative orientation of the magnetic 
cores. An expression for this energy is obtained using the 
magnetostatic field experienced by one of the dots due to 
the other. Details of these calculations are included in 
appendix. The resulting expression is 
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E^iS) = ±— I dk{l-e-^")Jo{kS) 



-ft 

j Mkp)'f 



,{p)pdp 
(10) 

Since the size of the core will be obtained by energy min- 
imization, in the previous equation we have considered 
that the two interacting dots exhibit the same magnetic 
profile, that is, their cores are identical. Equation pU]) 
allows us to write the interaction energy between two 
dots as E^{S) = ±E{S), depending on the relative ori- 
entation of the cores. Note that if the core size increases, 
i.e. niz grows, then E{S) increases. Also we always have 
E{S) > 0, and the interaction energy between two dots 
is always greater than zero {E~^{S) > 0) if the cores are 
oriented in the same direction, which causes a shrink- 
ing of the vortex core, in agreement with micromagnetic 
simulations [iOl . On the other hand, if the two cores are 
oriented in opposite directions, then the interaction term 
is negative {E^ (5) < 0) and the vortex core expands to 
lower the total energy. 

Substituting our expression for (p), Equation ([4|. 
we obtain 

i^(.) = 2-+Vr^(n + l)^ 

00 



Using equation (jlip and adding up contributions over 
the entire array we obtain the expression for the total 
interaction energy of the N x N square array as 

N-l 

Ei,{N) =2NY,{N- p){±lYE{pd) 



N-l N-l 



(12) 



This equation was previously obtained by Laroze et al, 
[1^ and has been used to investigate the magnetostatic 
coupling in arrays of magnetic nanowires. 

At this point we need to specify the value of the pa- 
rameter n for the core profile defined by equation ([?]). 
In a previous work, Landeros et al [l3| showed that the 
magnetic vortex core can be well described for almost 
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any value of n > 1. We choose n = 4, as explained in 
[14], and finally obtain the following expression for the 
self- energy, 



E.eif = 0.02987r6-^ 
with 



2.266J 
(13) 



1 



63a; 
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(256 + 384x2 ^ 57g^4 



+600a;^ + 350a:^ - 256^21 (-9/2, 1/2, 6, ~4:x'^)) , 



where F2i(a, 6, c, z) is a hypergeometric function. 
Using n = 4 in Eq. (Ilip we obtain 



E{s) = 294912 



^(1 
„io ^ 



(14) 

Finally, the total energy of the array is calculated as 
= N^Eseif+ E^^ , with E^eif given by equation ^ 



and E^f. given by equations and p^ . 

To determine the vortex core magnetization, we have 
to minimize E^ with respect to b. Note that the only 
term in the expression for the energy which depends on 
the radius r of the dot is Eseif (see equation ([T^ V How- 
ever, the derivative of Egg^if with respect to b is indepen- 
dent of r, leadin g to a core size that is independent of the 
dot radius, l2ll| This follows from the fact that the 
external region of the dot (a perfect vortex) does not in- 
teract with the core (apart from the exchange interaction 
across the interface between the two regions) . That is to 
say, the equation for b that minimizes the total energy of 
the vortex configuration is independent of r. 
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IV. RESULTS AND DISCUSSION 

We are now in position to investigate how the core ra- 
dius is affected by the interdot magnetostatic coupling. 
In order to obtain the value of 6, the energy E^ must be 
minimized for fixed h, d and N. Since the dipolar inter- 
action is long-ranged, an increment in N leads to an in- 
crease in the dipolar field felt by each dot and, therefore, 
to a change in the core radius until a certain asymptotic 
value is reached. Figure 3 illustrates this effect for arrays 
of dots with r — 6, d — 12.5 and 14, and h = 12 (figure 
3(a)) and h = 8 (figure 3(b)), in configurations P and 
AP. For both values of h we observe that for = 12 
the core radius is almost at its limiting size. From our 
results we observe an increase in the core radius with the 
number of dots in the AP configuration, and the oppo- 
site behavior for the P configuration, in agreement with 
micromagnetic simulations. [20[ For a given value of N, 
AP arrays with smaller values of d present larger core 
radius due to the preference for the AP ordering. The 
antiparallel alignment has the lowest dipolar energy, so 
the core radius increases with the number of dots in this 



FIG. 3: Effective core radius c^ff vs A'^^ for r — 6: (a) h = 12 
and (b) h = S. Open symbols correspond to an antiparalfef 
(AP) ordering and fulf symbofs to a paraffef (P) ordering. 
Circfes correspond to d = 12.5 and triangles to d = 14. 



configuration. On the other hand, in the P configuration 
the parallel coupling between the cores increases the in- 
teraction energy, and then the core region is reduced in 
order to decrease the interaction. 

The dipolar energy may also be varied if the interdot 
distance is changed. This effect is depicted in figure 4 
for the AP and P configurations. While in the antipar- 
allel arrays the size of the core rapidly reaches the value 
for isolated dots; in the parallel configuration, the effect 
of the interdot interaction is relevant for longer interdot 
distances. 

Now, the infiuence of the dot height is more subtle. 
Figure 5 shows the behavior of Ceff vs h for arrays of 
dots with r = 6, d = 12, with A^^ = 4 and 64, in the P 
and AP configurations. For the isolated dot (blue solid 
line) , the effective core follows approximately the relation 
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FIG. 5: Effective core radius c^jj vs h for dots with d — 12, 
r = 6. Open symbols correspond to an antiparallel (AP) 
ordering and full symbols to a parallel (P) ordering. Stars 
correspond to A''^ = 4 and circles to A'^^ = 64. 



function. 

The interaction energy between two dots with full mag- 
netization along their axis can be obtained from equation 
(14) using mz{p) — 1 and gives 



^(1 



(16) 



FIG. 4: Effective core radius Cg^j vs d for r = 6: (a) /i = 12 
and (b) h = S. Open symbols correspond to an antiparallel 
(AP) ordering and full symbols to a parallel (P) ordering. 



Ceff^ 1.228 + 0.291/i°"7. 

For an isolated dot, a transition from the vortex con- 
figuration to a complete ferromagnetic ordering along the 
dot axis is observed as h increases, as shown in the phase 
diagrams presented in [14]. From the point of view of 
the core radius, this transition may be seen as a slow 
and continuous increase in Ceff with h, until Ceff ~ r. As 
the dots interact in the array this behavior may change. 
To illustrate this point we calculate the transition line 
from the out-of-plane uniform state to the vortex-core 
state configuration in arrays of 4 and 64 dots. The self- 
energy for the out-of-plane uniform (u) state has been 
presented in [s^l and reads 



^self - 



1 + 



8r 



21 



4^2 
1^ 



(15) 



where -F21H = -F2i[— 1/2, 1/2, 2, is a hypergeometric 



as presented by Beleggia et al. [3l!| The total energy of an 
array with out-of-plane uniform magnetization is given 
by 



self ' 



E\ 



u± 



(17) 



with the above expressions for the self-energy (equation 
(fT5| l and the interaction energy {E"^^) given by equa- 
tion with the magnetostatic coupling of the two dots 
given by equation (|16p . The transition lines between the 
vortex-core and uniform out-of-plane magnetic states are 
obtained by equating equations ([7]) and p7|) and are rep- 
resented in figure 6. 

As we expected, the transition line shifts to lower radii 
(and greater heights) in the parallel array and to bigger 
radii (and lower heights) in the antiparallel array. This 
behavior can be understood by analyzing the interdot 
dipolar coupling between two dots, E{s). The absolute 
value of the interaction energy increases with the core 
size, so E^{s) > E{s), as the uniform out-of-plane (u) 
configuration can be seen as a vortex-core with an infinite 
core radius, and the magnetostatic coupling for uniform 
out-of-plane magnetization is stronger than the coupling 
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a transition from vortex to a full ferromagnetic state is 
hindered, while in the AP configuration the interaction 
favors the transition to a full AP ferromagnetic state. 
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APPENDIX 



FIG. 6: Dipolar induced shifts in the transition line that sepa- 
rates the vortex and the out-of-plane uniform magnetic states. 
We have fixed the ratio r/d = 0.48. Open symbols correspond 
to an antiparallel (AP) ordering and full symbols to a parallel 
(P) ordering. Stars correspond to A'^ = 2 and circles to A' = 8. 



of two vortex-cores. Therefore, the transition line in a 
parallel array shifts to the left, because the interactions 
in an array of vortex-core dots are less than the interac- 
tions in a uniformly magnetized P array. On the other 
hand, in the AP arrays, the interactions are negative and 
the uniformly magnetized AP configuration is favorable, 
so the complete ferromagnetic ordering is reached for a 
smaller value of h. 



V. CONCLUSIONS 

We have examined the influence of dipolar interactions 
in the effective core radius of dots in a vortex configura- 
tion, placed in an X A^ square array. Dipolar coupling 
among dots occurs via core interaction, so we consider 
two types of relative alignment between the cores: paral- 
lel, and antiparallel nearest neighbors. Whenever a cer- 
tain array configuration lowers the interaction energy, the 
core region expands, as occurs in AP interactions. The 
opposite behavior occurs for configurations increasing the 
interaction energy, as the P case. These two orderings 
represent a demagnetized configuration, the antiparallel 
case, and a saturated one, corresponding to the parallel 
case. Effects of the interaction between the dots in the 
core size have been investigated by varying the number 
of dots in the array, the distance between them, and their 
heights. In all cases, an increase of the dipolar interac- 
tion energy leads to a decrease of the core radius. When 
we increase the height of the dots, in the P configuration 



The dipolar interaction energy between two identical 
dots i and j in the vortex core configuration is given by 



where Mi is the magnetization of the i-th dot and Uj 
is the magnetostatic potential of the j-th dot. For the 
vortex core model defined by Eq. ([3])we have V • M = 0, 
so the magnetostatic potential reduces to [3] 



27r 



,ip(4>-4>')rlrh' 



P=-oo (J 

M,ip')Jp{kp')p'dp' J Jp{kpj)[e~^^""'-^ ^ e-^''-]dk . 



(A18) 

In this expression we have used the expansionist] 

(A19) 

Here Jp{z) are Bessel functions of first kind. The an- 
gular integration gives us J^^ e^P'^'^~'^'^d(j)' = 27re*P'^(5p,o, 
leading to 



Uj {Pj, p'Mz ip') Mkp')dp' 



Mkpj) 



dk. (A20) 



Since the potential has no dependence on (j)^ the expres- 
sion for the energy reduces to 

H 2n R 



{p)^^^^^g^pdpd^ dz , (A21) 
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and using Eq. (|A20p it is straightforward to obtain 

27r oo 





R 



dk I Ja{kpj)M^{p)pdp 
Jo(fcp')M.(p')(l - e-''")p'dp' . (A22) 



We need to evaluate the potential due to dot j on dot 
i a distance S apart. Then we have to relate the radial 
coordinates of both dots through the relation 



Pj - ^Jp-' + S^~2pScos{^ + P), 
with /3 an arbitrary angle. Using the following 



identityjH 



Jo(fcpj) = Jo (k^p"^ + 5*2 - 2pS cos{4> + /?)) 

OC 

= e^Pi^+'^)jp(kp)Jp{kS) 



in the expression for the energy (Eq. IA22[) . after the 
angular integration, we obtain 



E^[S]^±^ I dk{l-e~''")Jo{kS) 



LO 



jQ{kp)m^{p)pdp 
(A23) ■ 
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